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Abstract 

We discuss recent developments in the study of quantum wavefunctions and transport in classically ergodic sys- 
tems. Surprisingly, short-time classical dynamics leaves permanent imprints on long-time and stationary quantum 
behavior, which are absent from the long-time classical motion. These imprints can lead to quantum behavior 
on single-wavelength or single-channel scales which are very different from random matrix theory expectations. 
Robust and quantitative predictions are obtained using semiclassical methods. Applications to wavefunction in- 
tensity statistics and to resonances in open systems are discussed. 
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1. Introduction 

The structure of quantum wavefunctions and 
the closely related problem of quantum trans- 
port in classically non-integrable systems have re- 
ceived much attention recently from a variety of 
physics communities. Questions concerning the 
quantum behavior of systems with a generic clas- 
sical limit are of course of great fundamental in- 
terest; they are also very relevant not only for 
nanostructure and mesoscopics experiments [Q, 
but also for understanding phenomena in areas 
as diverse as atomic physics molecular and 
chemical physics microwave physics nu- 
clear physics and optics |^]. Combined with 
knowledge about spectral properties, wavefunc- 
tion information can be used to address conduc- 
tance curves, susceptibilities, resonance statistics, 
and delay times in ballistic quantum dots. Sim- 
ilar wavefunction and transport issues arise in 
other fields, in the study of resonance statistics 
in microwave cavities, photoionization cross sec- 
tions, chemical reaction rates, spectra of Rydberg 
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atoms, lifetimes and emission intensities for reso- 
nant optical cavities, and S'-matrix properties in 
many systems. 

In the classically integrable case, quan- 
tum wavefunctions are known to be asso- 
ciated with the invariant tori of the cor- 
responding classical dynamics, satisfying the 
Einstein-Brillouin-Keller (EBK) quantization 
conditions |^]. Classical-quantum correspon- 
dence in the ergodic case is, however, more sub- 
tle. Here, the typical classical trajectory uni- 
formly visits all of the energetically available 
phase space, so naively the typical quantum wave- 
function should also have uniform amplitude over 
an entire energy hypersurface, up to the in- 
evitable (Gaussian random) fluctuations. Such 
behavior follows directly from Random Matrix 
Theory (RMT), which has been proposed by Bo- 
higas, Giannoni, and Schmit Q to be the proper 
description of quantum chaotic behavior in the 
semiclassical limit (i.e. in the limit where the de 
Broglie wavelength A becomes small compared to 
the system size). A similar conjecture by Berry j^] 
states that a typical quantum wavefunction in 
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this same limit should look locally like a random 
superposition of plane waves of fixed energy, with 
momenta pointing in all possible directions. 

RMT @ , a natural quantum analogue of clas- 
sical ergodicity, turns out to describe well spectral 
properties on small energy scales (e.g. the distri- 
bution of nearest neighbor level spacings in the 
spectrum), but does not always provide a valid 
description of wavefunction structure and trans- 
port properties. Scars ||l^ provide one of the 
most visually striking examples of strong devi- 
ation from RMT wavefunction behavior; other 
examples include slow ergodic systems and 
Sinai-type systems In all these cases, non- 

RMT behavior can be quantified, semiclassically 
predicted, and observed. 

Section ^ addresses the problem of wavefunc- 
tion ergodicity in broad terms. We are interested 
both in general questions concerning the implica- 
tions of global classical properties (such as ergod- 
icity or mixing) on quantum behavior and also 
in the way in which specific classical structures 
(such as unstable periodic orbits) may leave their 
imprints on the quantum eigenstates. Several ex- 
amples, including that of Sinai-type systems, are 
discussed in Section 2^. Sections § and ^ focus 
on the scar phenomenon and quantitative predic- 
tions. The overview format of this presentation 
requires us to omit most derivations; references 
to more detailed discussions in the literature can 
be found throughout. 

2. Ergodic wavefunction structure and 
quantum transport 

2.1. Coarse-grained vs. microscopic ergod- 
icity 

We must distinguish between two ways of ex- 
tending classical notions of ergodicity to the 
quantum case fl^ . First, we may consider wave- 
function intensity integrated over a classically de- 
fined region 7^, in the regime where the wave- 
length A becomes small compared to the size of TZ. 
As can be shown rigorously [Q , in this limit the 
integrated intensity approaches a constant (equal 
to the area of 7?. as a fraction of total phase 
space) for almost all wavefunctions. The result 
requires only long-time classical ergodicity and 



quantum-classical correspondence at short times; 
a half-page physicists' derivation can be found in 
|l5|| . We note that this macroscopic or "coarse- 
grained" type of quantum ergodicity is clearly 
implied by RMT, but is in fact a much weaker 
condition. Quantum wavefunction structure on 
coarse-grained scales is relevant for studying con- 
ductance and conductance fluctuations through 
wide, multichannel leads, for fast decay processes, 
and generally for analyzing open systems in the 
overlapping resonance regime. 

RMT, on the other hand, is a prediction about 
wavefunction uniformity at the quantum scale, 
i.e. on the scale of a single wavelength (or sin- 
gle momentum channel, or most generally at the 
scale of a single ?i-sized cell in phase space). This 
kind of uniformity is a much stronger condition 
than coarse-grained ergodicity, and several exam- 
ples will be given later in this section where mi- 
croscopic quantum ergodicity is violated in clas- 
sically ergodic systems. Wavefunction structure 
at the quantum scale is relevant, obviously, for 
transport through narrow (or tunneling) leads. 

2.2. Measures of microscopic ergodicity 

First we must define quantitative measures 
of ergodicity or localization at the microscopic 
scale Let \n) be an iV— dimensional basis of 
eigenstates, and \a) some localized test basis, e.g. 
position, momentum, or phase space Gaussians, 
as is physically appropriate in a given system. For 
example, in discussing Anderson-type lattice lo- 
calization, we may choose our test basis |a) to 
be the position basis, whereas in scattering prob- 
lems plane waves may be a more natural choice. 
Sometimes we take \a) to be the eigenstates of a 
zeroth-order Hamiltonian Hq or of a zeroth-order 
scattering matrix 5*0. We are then interested in 
the wavefunction intensities Pan = and 
their correlations. For simplicity of presentation 
we assume there are no conserved quantities, so 
classically nothing prevents each eigenstate \n) 
for overlapping equally with each \a) (of course, 
the formalism generalizes naturally to the more 



general case of classical symmetries, see |16 ). We 
adopt the normalization convention 



{Pan) = 1: 



(1) 
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where the averaging is done over test states |a), 
over wavefunctions |n), or over an appropriate en- 
semble of systems. 

The first nontrivial moment of the Pan distri- 
bution is the inverse participation ratio: 

1 ^ 

I™° = ]^E^->1, (2) 

71=1 

which measures the mean squared wavefunction 
intensity at \a), or, alternatively, the inverse frac- 
tion of wavefunctions having significant intensity 
at \a). IPR„, the mean squared intensity of a sin- 
gle wavefunction \n) averaged over position |a), is 
defined analogously: it measures the inverse frac- 
tion of phase space covered by |n). A global IPR 
measure may also be conveniently defined: 

N N 

lPR=-J2lPRa^JjY.^PRn>l, (3) 

a— 1 n— 1 

and provides a simple one-number measure of the 
degree of localization in a quantum system. An 
IPR of unity corresponds to perfect ergodicity 
(each wavefunction having equal overlaps with all 
test states), while IPR = N indicates the greatest 
possible degree of localization (the wavefunctions 
|n) being identical with the test states \a)). 

Besides being the first nontrivial moment of the 
intensity distribution, the IPR measure is useful 
because of its connection with dynamics. For ex- 
ample, IPRq is proportional to the averaged long- 
time return probability for a particle launched ini- 
tially in state |a): 

IPR, = 7V hm i^'|(a|e-^*/^|a)p. (4) 

t=0 

It is intuitively clear that enhanced long-time re- 
turn probability is associated with increased lo- 
calization. We may generalize IPRq to a trans- 
port or wavefunction correlation measure between 
local states \a) and \b): 

P^, = f: = lun^ ^ j:\{a\e-^^^^^\b)\' .(5) 

ji=l t=o 

Of course, the mean J^t Pab = 1 by unitarity; 
the simplest nontrivial measure of transport effi- 



ciency is thus 

Qa = ^Y.Pa,>^, (6) 

b 

which measures the inverse of the fraction of 
phase space accessible from \a) at long times. 
Qa = l for all \a) (the RMT result for oo) 
indicates perfect long-time transport, and vanish- 
ing wavefunction correlations. 

2.3. Examples 

The Random Matrix Theory description is free 
of all dynamical information about the system un- 
der study and thus can hardly be expected to pro- 
vide a correct statistical description of all quan- 
tum behavior. It serves, however, as a very useful 
baseline with which real quantum chaotic behav- 
ior may be compared. In RMT, the wavefunction 
intensities Pan are squares of Gaussian random 
variables and so are drawn from a distribution 
(of one degree of freedom for real overlaps {a\n) or 
two for complex overlaps), with a mean value of 
unity. This easily leads to IPRrmt = 3 in the real 
case or 2 in the complex case. Notice that even in 
RMT, quantum fluctuations cause wavefunctions 
to be less ergodic than the classical expectation 
IPRcias = 1- By quantum localization, however, 
we always mean fluctuation in the intensities in 
excess of what would be expected from a Gaus- 
sian random model, i.e. IPR > IPRrmt- Ex- 
amples abound of such anomalous quantum be- 
havior in classically ergodic systems, and several 
are described below. We also mention here that 
in RMT, the channel-to-channel transport effi- 
ciency Pab approaches unity for all channels |a), 
1 6) in the N ~> oo semiclassical limit. 

(i) Scarring is the anomalous enhancement or 
suppression of quantum wavefunction intensity on 
the unstable periodic orbits of the correspond- 
ing classical system. This localization behavior 
is perhaps surprising from a naive classical point 
of view, since in the time domain (see Eq. ^ it 
implies an enhanced long-time return probability 
for a wavepacket launched on an unstable peri- 
odic orbit. Paradoxically, this is in contrast with 
the classical behavior, where a probability dis- 
tribution spreads itself evenly over the entire er- 
godic space at long times and retains no memory 
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of its initial state. Quantum long-time dynamics, 
which of course contains phase information, thus 
retains a much better memory of the short-time 
classical behavior than does the long-time clas- 
sical dynamics, for arbitrarily small values of h. 
Specifically, for a wavepacket \a) optimally ori- 
ented on an orbit of instability exponent /?, we 
have inll 



IPRsc 



E 



m— — oo 



cosh /3m 



IPRrmt (7) 



/3 



IPR.RMT ^ IPR.RMT , 



(8) 



where the limiting form is valid for weakly un- 
stable orbits (/3 ^ 1). A finite fraction of order 
(3 of all wavefunctions arc scarred: they have in- 
tensity 0{f3~^) on the periodic orbit compared 
with the mean intensity, while most of the re- 
maining wavefunctions are antiscarred, and have 
intensities on the orbit as small as exp {—7t^/2(3) 
compared with the mean. The source of this lo- 
calization effect will be outlined in Section ^. The 
transport measure Pab is also affected by scarring: 
if both leads \a) and \b) are located near periodic 
orbits, Pab can be enhanced or suppressed de- 
pending on whether the two orbits are "in-phase" 
or "out-of-phase" in the range of energies con- 
sidered. These long-time transport results are 
?i— independent and are obtained directly from 
the short-time linearized classical dynamics. 

(ii) The tilted wall billiard (a rectangular box 
with one wall tilted relative to the other three), 
is a classically ergodic system in which the wave- 
functions (as measured by the IPR) become less 
and less ergodic in the classical h (or A ^ 0) 
limit . (This kind of behavior is possible be- 
cause the h ~^ limit fails to commute with 
the infinite time i — > oo limit.) Specifically, for 
wavefunctions expressed in channel (momentum) 
space, we have 



IPRxiitcd - A-i/VlogA-i ^ oo. 



(9) 



This anomalous wavefunction behavior at the 
single- channel scale is entirely consistent with er- 
godicity on coarse-grained scales; coarse-graining 
over ~ A^^/logA^^ channels is required to re- 
trieve the classical behavior. Quantum transport 



in these systems is also anomalous and is domi- 
nated by diffractive effects (due to the slowness 
of classical phase space exploration). 

(iii) As a final example we mention the Sinai 
billiard, a paradigm of classical chaos jl^. This 
consists of a circular obstruction of diameter d ^ 
A placed inside a rectangle. Again using momen- 
tum channels for the test basis |a), we obtain 



IPRsi, 



(logA-i)/d 



(10) 



for fixed d in the A ^ limit. Again, we see 
nonergodic A — > wavefunctions in a classically 
ergodic system. The distributions of wavefunc- 
tion intensities Pan and IPR's (IPRa and IPRn, 
see Eq. |^) all display power-law tails (in contrast 
with the exponential tail prediction of RMT). 
Transport is similarly anomalous: for the typi- 
cal channel \a), Qa ^ d^^ (see Eq. ^); i.e. at long 
times a given channel is coupled to only a fraction 
~ d ^ 1 of all other channels. 

3. Scars of quantum chaotic wavefunctions 

3.1. Phenomenology and conceptual issues 

Scars have been observed experimentally and 
numerically in a wide variety of systems. These 
include semiconductor heterostructures, where 
scars are observed to affect tunneling rates from 
a 2D electron gas into a quantum well mi- 
crowave cavities p9[ |, the hydrogen atom in a 
uniform magnetic field |2^], acoustic radiation 
from membranes pl| , and the stadium billiard, 
for which wavefunctions as high as the millionth 
state in the spectrum can be numerically ana- 
lyzed |2^. Nevertheless, confusion over the def- 
inition and measures of scars (and the dearth of 
quantitative predictions) have until recently led 
some to question the existence of a scar theory, 
and even to doubt the survival of the phenomenon 
in the classical limit. Recent theoretical develop- 
ments jl^] enable us to make robust, quantita- 
tive predictions about how strongly a given orbit 
will be scarred and how often, as a function of 
energy and other system parameters; real com- 
parison with numerical and experimental data is 
therefore made possible. We also note that statis- 
tical and quantitative data, rather than anecdotal 
(in the form of wavefunction plots) is necessary as 
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a test of scar theory or of any other theory of a 
locahzation phenomenon, since considerable fluc- 
tuations of wavefunction intensity occur even in 
the context of R MT (see discussion in first para- 
graph of Section 2.3). 

We address some common misconceptions and 
summarize key facts about the scar phenomenon: 
(i) Scarring is associated with unstable periodic 
orbits, not with stable or marginally stable ones, 
though these of course also do attract wavefunc- 
tion intensity; a qualitative difference arises be- 
cause of classical-quantum noncorrespondence in 
the unstable case (Section 2.3(i)). (ii) Weak scars 
are not always visible to the naked eye; scars are 
defined and measured according to a statistical 
definition, (iii) Scarring predictions are robust, 
valid even when the exact dynamics is not known 
well enough to allow individual eigenstates to be 
determined either quantum mechanically or semi- 
classically. (iv) The amount of scarring associated 
with individual eigenfunctions varies significantly 
from state to state, but in accordance with a the- 
oretically predicted distribution, (v) If an opti- 
mally chosen phase space basis is used, the typ- 
ical intensity enhancement factor as well as the 
full distribution of scar intensities can be given 
as a function of the instability exponent /3. 

3.2. Short-time effects on stationary prop- 
erties 

We define the autocorrelation function for test 
state la): 



A{t) = {a\ exp[-iHt] 



(11) 



the Fourier transform of which is the local density 
of states (LDOS) 



SiE) = PanS{E - En) . 



(12) 



Then, if we know the statistical properties of the 
return amplitude A{t), we also know the statis- 
tical properties of the LDOS at |a), and specifi- 
cally of the wavefunction intensities Pan (see e.g. 
Eq. Q). Furthermore, given information about 
the return amplitude A{t) for short times only 
(say for \t\ < Tq), we immediately obtain the 
LDOS envelope S'smooth (^') , which is the energy- 
smoothed (on scale Tl/Tq) version of the true line- 



spectrum S{E). So large short-time recurrences 
in A{t) get 'burned into' the spectrum. 

Specifically, if \a) is a minimum uncertainty 
wavepacket optimally placed on a periodic orbit 
of period P, then at integer multiples t — mP 
of this period, by Gaussian integration we easily 
obtain ITlll 



^short(^) — ex\i[imS /%]/ \/ cosh /3m , 



(13) 



where the classical orbit has action S and insta- 
bility exponent /?. For small /3, the Fourier trans- 
formed spectral envelope 5'smooth(-E) has bumps 
of width scaling as (3/P and height scaling as /3~^ . 
Now after the mixing time (^ /3~^logA^), es- 
caping probability starts returning to the origin, 
leading to fluctuations in the spectral envelopes, 
and eventually to discrete delta-function peaks 
(Eq. |l2|). The nonlinear scar theory allows the 
long-time returning amplitude to be analyzed as 
a homoclinic orbit sum [ pl| , and leads to the pre- 
diction that the spectral intensities are given by 
variables multiplying the short-time envelope: 



Pa 



(14) 



where the ran follow a Gaussian random distri- 
bution of variance unity. 

4. Quantitative predictions for scar effects 

4.1. Wavefunction intensity statistics 

As a simple application, we may compute the 
distribution of wavefunction intensities on a pe- 
riodic orbit of instability exponent (3 For 
complex wavefunctions, the probability to have 
an intensity exceeding the mean by a factor x is 
given by exp(— x) in RMT; on a periodic orbit the 
tail of the distribution is instead given by 



P{P,x) = C(3{l3x) 



(15) 



where G and Q are known numerical constants. 
Similarly, the tail of the overall intensity distri- 
bution (sampled over the entire phase space) is 
given by 



P{p,x) =C'/3h{Px)-^/^e-^^^^ , 



(16) 



where /? is now the exponent of the least unstable 
periodic orbit. These two results are illustrated 
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Figure 1. Cumulative wavefunction intensity dis- 
tribution (a) as measured by a Gaussian test state 
optimally centered on a periodic orbit with insta- 
bility exponent /? = log 2, plotted as the upper 
thick curve with scar theory prediction for the 
tail (Eq. |l|) given by dashed curve, and (b) av- 
eraged over the entire phase space of size 200ft., 
plotted as lower thick curve with theory for the 
tail (Eq. |l6|) given by solid curve. The dotted line 
is the RMT exponential prediction (after ||2^). 
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Figure 2. The long-time enhancement of the 
probability to remain in a weakly open system 
when the tunneling lead is placed on a periodic 
orbit is plotted as a function of the instability ex- 
ponent /3 of the orbit. Data is shown for system 
sizes = 120 (plusses) and N — 240 (squares). 
The N ^ oo theoretical prediction (Section ^^ ), 
growing exponentially with decreasing /3, is shown 
as a solid curve. For large (3, the enhancement 
factor converges to 1, the RMT prediction (af- 
ter [pi). 



in Fig. |l|. Finally, for an ensemble of systems, a 
power-law tail is obtained, in contrast with the 
exponential prediction of RMT, and large inten- 
sities have been observed numerically with fre- 
quency exceeding by ~ 10'^'' the RMT predic- 
tions. 

4.2. Enhancement in probability to remain 

The above methods can also be applied to the 
open systems case [Q. Consider the probability 
to remain inside a chaotic quantum well coupled 
to the outside via a tunneling lead. For a lead 
optimally placed with respect to a periodic orbit 
of exponent /3 1, the long-time probability to 
remain in enhanced by a factor exp(7r^/2/3) (see 
Fig. H), compared with the RMT expectation. 
(The en han cement is due the antiscarred states 
(Section 2.3(i)), which have very small coupling to 
the lead.) Of course, the classical probability to 
remain is exponential and independent of lead po- 
sition for a chaotic system with a narrow lead. So 
once again we see long-time quantum mechanics 
retaining an imprint of short-time classical struc- 
tures which is absent from long-time classical be- 



havior. The analysis can of course be extended 
to the study of conductance peak statistics in two 
lead systems, where one or both leads are located 
near short periodic orbits (compare with the dis- 
cussion of Pab statistics in Section |2^(i)). 

5. Conclusions 

We have seen that wavefunction structure and 
transport in classically ergodic systems, including 
the paradigmatic Sinai billiard system, can dif- 
fer greatly from RMT expectations, and in fact 
can deviate further and further from RMT in the 
semiclassical limit. This non-ergodic quantum 
behavior at the scale of single wavelengths or sin- 
gle quantum channels can be quantitatively and 
robustly predicted using only short-time classical 
information. This anomalous small-scale quan- 
tum behavior is also entirely consistent with er- 
godicy on coarse-grained scales, as studied by 
Schnirelman, Zelditch, and Colin de Verdiere. 

Scarring is a fascinating example of the influ- 
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ence of identifiable classical structures on sta- 
tionary quantum properties (e.g. eigenstates) 
and long-time quantum transport in classically 
chaotic systems. Short unstable periodic orbits 
leave a strong imprint on the long-time proper- 
ties of the quantum chaotic system, even though 
classical dynamics loses all memory of these struc- 
tures at long times. Scar theory makes robust and 
quantitatively verified predictions about proper- 
ties such as the wavefunction intensity distribu- 
tion in a chaotic system, including the power- 
law tail observed after ensemble averaging. A 
lead centered on an unstable periodic orbit has 
been shown to produce many exponentially nar- 
row quantum resonances, even though the decay 
times are very long compared with all other time 
scales in the problem. Enhancement factors of 10 
to 100 in the total probability to remain are easily 
observed for moderate values (1.0 to 0.5) of the 
instability exponent. 
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